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Abstract —Constant envelope (CE) precoding is an appealing 
transmission technique which enables the realization of high 
power amplifier (PA) efficiency. For CE precoding in a single- 
user multiple-input single-output (MISO) channel, a desired 
constellation is feasible at the receiver if and only if it can be 
scaled to lie in an annulus, whose boundaries are characterized 
by the instantaneous channel realization. Therefore, if a fixed 
receiver constellation is used for CE precoding in a fading 
channel, where the annulus is time-varying, there is in general 
a non-zero probability of encountering a channel that makes 
CE precoding infeasible, thereby causing a high probability of 
error. To tackle this problem, this paper studies the adaptive 
receiver constellation design for CE precoding in a single-user 
MISO fiat-fading channel with an arbitrary number of antennas 
at the transmitter. We first investigate the fixed-rate adaptive 
receiver constellation design to minimize the symbol error rate 
(SER). Specifically, an efficient algorithm is proposed to find the 
optimal amplitude-and-phase shift keying (APSK) constellation 
with two rings that is both feasible and of the maximum 
minimum Euclidean distance (MED), for any given constellation 
size and instantaneous channel realization. Numerical results 
show that by using the optimized fixed-rate adaptive receiver 
constellation, our proposed scheme achieves significantly im¬ 
proved SER performance over CE precoding with a fixed receiver 
constellation. Furthermore, based on the family of optimal fixed- 
rate adaptive two-ring APSK constellation sets, a variable-rate 
CE transmission scheme is proposed and numerically examined. 

Index Terms —Constant envelope (CE) precoding, adaptive 
receiver constellation, amplitude-and-phase shift keying (APSK), 
minimum Euclidean distance (MED). 

1. Introduction 

Improving the power efficiency of radio frequency (RF) 
power amplifiers (PAs) reduces the energy consumption of 
wireless communication systems. From the perspective of 
power efficiency maximization, the most favorable input sig¬ 
nals for PAs are constant envelope (CE) signals. Specifically, 
since CE input signals have the lowest possible peak-to- 
average power ratio (PAPR), the backoff required for PA 
operation is minimized, and hence the power efficiency is max¬ 
imized 1^. It is also worth noting that CE input signals impose 
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lower requirement on the dynamic range of the PAs than their 
non-CE counterparts, thereby requiring less expensive PAs to 
be used in practice. 

To realize CE input signals, the complex baseband signal at 
each transmit antenna is required to have constant amplitude, 
and information can only be modulated in the transmitted 
signal phase. This stringent constraint, however, induces new 
signal processing challenges. For single-antenna channels, 
various CE modulation techniques have been studied, e.g., 
continuous phase modulation (CPM) ||2L CE-OEDM (orthog¬ 
onal frequency division multiplexing) |[^, etc. On the other 
hand, for multi-antenna channels, CE precoding was recently 
proposed and investigated in 0-191. Specifically, for the case 
of single-user multiple-input single-output (MISO) channels, it 
has been shown in 0 0 that by varying the signal phases at 
different transmit antennas, the noise-free signal at the single¬ 
antenna receiver always lies in an annular region (between 
two concentric rings), whose boundaries are determined by 
the instantaneous channel realization and per-antenna transmit 
power. In addition, low-complexity algorithms were proposed 
in 0 for CE precoding to achieve a nonlinear mapping 
from any desired received signal point within the annulus 
to the corresponding transmitted signal phases based on the 
instantaneous channel state information (CSI)[^ Moreover, 
for multi-user large-scale MISO downlink systems ef¬ 
ficient CE precoding algorithms were developed in 0, 0 
for frequency-fiat channels and were shown to guarantee an 
arbitrarily low multi-user interference (MUI) power at each 
user receiver with a sufficiently large number of transmit 
antennas. The work in 0,0 has been extended to the case 
with frequency-selective channels in 0, where the transmitted 
signal phases for consecutive channel uses are jointly de¬ 
signed. Eurthermore, 0 extended the CE precoding scheme 
proposed in 0 to the case with an additional restriction on the 
phase difference between consecutive channel uses, in order to 
prevent the potential spectral regrowth caused by abrupt phase 
changes. Although the per-antenna CE constraint is clearly 
more restrictive than conventional average-based sum power 
constraint (SPC) and per-antenna power constraint (PAPC) 
(see, e.g., (Tg-(T7)) , it was shown that with M transmit 
antennas, an array power gain of 0{M) is still achievable 
with the CE precoding schemes proposed in 0 - 0 - 

^It is worth noting that there is another line of research on prototyping and 
performance analysis of the so-called outphasing or LING (linear amplification 
with nonlinear components) technique (see, e.g., (^, ( 0 ), which can be 
viewed as a special case of CE precoding with two transmit antennas and one 
receive antenna. 
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However, note that even for the simplest case of the single- 
user MISO channel, a desired constellation at the receiver is 
feasible with CE precoding if and only if it can be scaled 
to lie in the annular region, i.e., all the signal points in the 
scaled constellation can be mapped back to CE signals at the 
transmitter. Therefore, for a fading channel where the annulus 
is time-var\dng, a fixed receiver constellation may not always 
be feasible]^ This can lead to severe performance degradation, 
which thus motivates this work on designing adaptive receiver 
constellations based on the instantaneous CSI. There are 
generally two approaches for adaptive receiver constellation 
design depending on whether the transmission rate or 
constellation size N is fixed or adjustable for the given 
application: 

• Eor delay-constrained systems that require fixed-rate trans¬ 
mission (e.g., real-time voice and video), N is fixed and the 
constellation S is adapted to channel conditions to minimize 
the receiver symbol error rate (SER); 

• Eor delay-tolerant systems that allow for variable-rate trans¬ 
mission, both N and S can be jointly adapted based on the 
CSI to maximize the average transmission rate subject to a 
given SER requirement at the receiver. 

To the best of our knowledge, neither approach in the above 
has been addressed in the literature for CE precoding. 

In this paper, we study the adaptive receiver constellation 
design for CE precoding in a single-user MISO fiat-fading 
channel with arbitrary number of antennas at the transmitter. 
Both cases of fixed-rate and variable-rate transmissions are 
considered. Our main contributions are summarized as follows: 

• Eirst, we derive the fixed-rate adaptive receiver constellation 
design that minimizes the SER at the receiver. This problem 
belongs to the class of circle packing problems that are 
known to be NP-hard fT^ . Hence, we approximate the exact 
SER by its union bound and assume that the constellation 
takes the form of amplitude-and-phase shift keying (APSK) 
so as to find a tractable solution. An efficient algorithm 
is proposed to obtain the optimal APSK constellation with 
two rings that is both feasible and maximizes the minimum 
Euclidean distance (MED) for SER minimization, given any 
constellation size and instantaneous CSI. Note that unlike 
existing works on APSK constellation design (see, e.g., p0| , 
1^ ) which consider an average power constraint of the 
signal points, our design is subjected to peak and minimum 
power constraints to ensure the constellation feasibility. It 
is also worth noting that for a given constellation size, our 
proposed design depends only on the ratio of the inner and 
outer radii of the annulus resulting from the instantaneous 
channel realization, and yields only a finite number of 
constellations, each corresponding to a certain range of 
this ratio. Therefore, the adaptive constellation set for any 
desired constellation size can be designed offline and stored 

^However, there are certain cases in which the CE precoding is always 
feasible regardless of channel fading. For example, PSK (phase shift keying) 
constellation is always feasible provided that the outer radius of the annulus 
is not zero. As another example, for a large-scale MISO channel with inde¬ 
pendent and identically distributed (i.i.d.) Rayleigh fading, the annulus was 
shown to become a disk region 0,0. therefore any receiver constellation 
is feasible. 


at the transmitter/receiver for real-time transmission. To 
further reduce the amount of memory required for storage, 
which increases with the number of constellations resulting 
from optimization, we also propose a suboptimal design of 
the adaptive two-ring APSK constellation set which consists 
of fewer constellations than the optimal set. 

• Next, for the case of fixed-rate transmission, we present 
numerical results which show that our proposed scheme 
significantly outperforms CE precoding with fixed receiver 
constellation in terms of average SER, which is due to our 
constellation adaptation and optimization tailored for CE 
precoding 1^ It is also shown by numerical results that for 
our proposed scheme, replacing the optimal constellation 
set with the suboptimal one only leads to small performance 
loss that is almost negligible. Moreover, the effect of imper¬ 
fect CSI at the transmitter (CSIT) on the performance of our 
proposed schemes is also investigated. 

• Einally, based on the family of optimal fixed-rate adaptive 
two-ring APSK constellation sets, we propose a variable- 
rate transmission scheme for CE precoding. According to 
the instantaneous CSI, the constellation size is selected as 
the maximum value from its feasible set that yields the 
SER union bound lower than a target value, to maximize 
the transmission rate. The performance of the proposed 
scheme is numerically examined in terms of average spectral 
efficiency and compared with that of another variable-rate 
CE transmission scheme based on the family of rectangular 
QAM (quadrature amplitude modulation) constellations. 


The remainder of this paper is organized as follows. Section 
presents the system model for CE precoding. Section 
ni| presents the problem formulation of fixed-rate adaptive 
receiver constellation design, then Section and Section V 
present the optimal solution and a suboptimal solution, respec¬ 
tively. Numerical results are shown in Section |V^ where the 
variable-rate transmission scheme is also presented. Einally, 
Section VII concludes the paper. 

Notations’. Scalars and vectors are denoted by lower-case 
letters and boldface lower-case letters, respectively, denotes 
the transpose of a vector z. ||z||i and ||z||oo are the /i-norm and 
loQ-norm of a vector z, respectively. | 2 ;| is the absolute value 
of a scalar 2 ;. denotes the space of m x n complex 

matrices. The distribution of a circularly symmetric complex 
Gaussian (CSCG) random variable with mean p and variance 
is denoted by and ~ stands for “distributed 

as”. Prob(-) denotes the probability. max{x,^} denotes the 
maximum between two real numbers x and y. \^\ denotes the 
cardinality of a set A'. A' U 3^ denotes the union of two sets 
A’ and 3^. 


H. System Model 

We consider a single-user MISO fiat-fading channel with 
M antennas at the transmitter and one antenna at the receiver. 
The baseband received signal is given by 

y = lV'x + n, (1) 

^For example, the proposed optimal 16-APSK with two rings has an MED 
gain over rectangular 16-QAM (0.5411 versus 0.4714) when the latter is 
feasible. 
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where x denotes the M x 1 transmitted signal vector; h G 
C^xi denotes the channel vector, which is assumed to be 
perfectly known at both the transmitter and receiver unless 
specified otherwise; n ^ CA/'(0,cr^) denotes the CSCG noise 
at the receiver. Assume a total transmit power denoted by P, 
which is equally allocated to the M antennas. Under the per- 
antenna CE constraint, the transmitted signal at each antenna 
is expressed as 


i.i.d. zero-mean CSCG random variables (i.e., i.i.d. Rayleigh 
fading channels), it can be shown that Prob( R > 3 ) “ 0-40 
i.e., rectangular 16-QAM constellation is infeasible with CE 
precoding for 40% of channel realizations. This will result 
in severe performance degradation, which thus motivates our 
design of adaptive receiver constellation for CE transmission 
based on the instantaneous CSI. 



where information is modulated in the transmitted signal 
phases Oi G [0,27r),i = 1,...,M. By varying 6>i’s, the noise- 

free received signal d = h^x = YliLi is shown 

in 0 and 0 to always lie in the following region: 




V= {deC:r <\d\<R}, (3) 

where R = .y/^||h||i and r = .y/^max{2||h||oo - ||h||i, 0}. 

As a result, any desired received signal constellation S (e.g., 
QAM) is feasible if and only if there exists a scaling factor 
a > 0 such that aS C V. In other words, S is feasible 
if and only if ^ < min|s|/max|s|. Eor any feasible S, the 

transmitted phase set {Oi} corresponding to any signal point 
s G 5 at the receiver and any a such that d = as can be 
readily obtained by the algorithms proposed in ||4||and 0> 
which are thus omitted for brevity. The channel inQ is thus 
equivalently represented by 

y = as n. (4) 

Eor illustration, Eig. shows the system diagram for CE 
precoding with M = 2 transmit antennas. Note that in order 
to maximize the received signal power and yet meet the 
feasibility condition, we should set a = R for any feasible 
constellation S with max|s| = 1, such that the signal point 

with the largest amplitude in S lies on the outer boundary of 
V at the receiver. 



Fig. 1. System diagram for CE precoding with M = 2. 

Eor MISO fading channel with time-varying h^s and as a re¬ 
sult both time-varying r and P, a fixed receiver constellation S 
may not always be feasible. Eor example, if 5 is a rectangular 
16-QAM constellation (i.e., minis I/maxi 51 = 4 ), it is feasible 

s s G*S 

only when ^ holds; otherwise, it is infeasible, as shown 
in Eig.|^ In particular, for the case where M = 2 and h^s are 


(a) Feasible case with ^ \ (b) Infeasible case with ^> \ 

Fig. 2. Feasibility of rectangular 16-QAM with different channel realizations. 


III. Problem Eormulation 


In this section, we formulate the problem of fixed-rate 
adaptive receiver constellation design for CE precoding. Our 
aim is to optimize the A^-ary constellation S with fixed N 
based on the instantaneous CSI to minimize the receiver SER. 
Note that for given N, our constellation design depends only 
on the ratio ^ resulting from /i^’s, but not on the exact 
values of or the value of M. Therefore, our proposed 
constellation design can be implemented offline and then used 
for real-time transmission. Also note that as will be shown later 
in the paper, there is only a finite set of constellations in our 
design for any given N, each corresponding to a certain range 
of the value Hence, the amount of memory required for 
storing the designed constellations at the transmitter/receiver 
is manageable. 

However, this constellation optimization problem is in gen¬ 
eral challenging, since it belongs to the class of circle packing 
problems that are known to be NP-hard m- Therefore, 
we make the following simplifications to obtain a tractable 
solution: 


1) Instead of considering the exact SER, P^, at the receiver, 
we minimize the union bound of SER given by P^ < 


(AT - 


12 


where we assume that S 


2 ) 


is an equiprobable signal set with max|s| = 1, and the 

sES 


maximum likelihood (ML) detection is used at the receiver 
to recover the signal point in S. This can be achieved by 
maximizing dmin, which denotes the MED between any 
two constellation points in S\ 

We assume that the constellation S follows an APSK 
structure, where signal points are distributed on concen¬ 
tric rings. APSK constellations are favorable due to the 
following reasons: first, they can be easily fitted into any 
annular region at the receiver by adjusting the radii of their 


^In this case, the cumulative distribution function (CDF) of ^ can be shown 
to be given by Prob(^ < x) = , x G [0,1]. 
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rings; moreover, since only a limited number of parame¬ 
ters suffice to fully characterize an APSK constellation, 
adaptive APSK can be designed and implemented with low 
complexity, as will be shown later in this paper. 


Consider an A^-ary APSK constellation S modeled by 

5= l = (5) 

where L > 1 denotes the number of concentric rings; pi 
denotes the radius of the Ith ring; Ni > 1 denotes the number 
of signal points that are uniformly spaced on the Ith ring, with 
Ni = N\uji denotes the phase offset of the Ith ring with 
respect to the reference phase cci = 0. We assume the L rings 
are indexed such that pL < Pl-i < < Pi- Thus, we have 

Pi = 1 since max|s| = 1 by assumption. 

seS 

Our objective is to maximize the MED of S by jointly 
optimizing L, {pi}iL 2 subject to all 

the constellation points lying in an annulus with inner and 
outer boundary radii given by ^ and 1, respectively. Let 
jC = {1,2,...,!/} denote the set of rings. We formulate the 
following optimization problem with given N and ^ as 

(PI) max dmin 

{Pz}z = 2’{^Tz = 2’^min 

S.t. dmin <du Ml e C (6) 

G^min < di^h, Ml,he jC, l^h (7) 
— < pL < pL-1 < • • • < P2 < 1 (8) 


L 


J2Ni=N 

1 = 1 

(9) 

Ni e Z+, yieC 

(10) 

uil G [0,27r), V/ G /^\{1} 

(11) 

L G Z+, 

(12) 


where is the set of positive integers; di denotes the intra¬ 
ring MED between any two signal points on the Ith ring; di^h 
denotes the inter-ring MED between any two signal points that 
are located on the Ith and hth rings, respectively. Erom ([^, di 
can be obtained as 



Note that when Ni = 1, di = oo and the corresponding 
constraint in ([^ of (PI) can be removed for ring 1. Using 
the cosine rule, we can also obtain 

di,h = \/pf + pl-‘^PiPhCi,h{Ni,Nh,uJi,uJh), I h, (14) 
where 

Cih{Ni,Nh,i0uL0h) = max cos 

m,n 

me{0,l,-,Nh-l}, ne{0,l,-,Ni-l}. (15) 


/ 27rn 27rm\ 


In this paper, we solve (PI) for the special case of L = 2|^ 
while our proposed method can also be extended to the more 
general case of L > 2 in principle. With L = 2, Problem (PI) 
is simplified as 


(P2) max dmin 

P2,<^2,N2,dmiri 


s.t. dmin - vGb} < 0 

(16) 

dmin - P2\/^2 < 0 

(17) 

G^min ~ a/i + P 2 “ 2p2C'i,2(A^2,^2) < 0 


(18) 

}<«<! 

(19) 

UI 2 € [0, 27r) 

(20) 

Ar2G{i,...,7V-i}, 

(21) 

f 1 - cos > 2 

with Ni = N - N 2 , Bi = 1 - 

for 


oc, Ni = 1 


I e {1,2}, and Ci, 2 (Al 2 , ^ 2 ) = max cos + W 2 - , 

m e (0,1 ,Ni - 1), n G (0,1, N 2 - 1). 

IV. Optimal Solution to Problem (P2) 

Note that Problem (P2) is a non-convex optimization prob¬ 
lem since the constraint in •El is non-convex due to the 
coupled p 2 and UJ 2 , and the fact that N 2 is an integer variable. 
In this section, we propose an efficient algorithm for solving 
Problem (P2). Specifically, we first optimize UJ 2 and p 2 for 
each given N 2 G {1,2, ...,A^ — 1}, and denote the correspond¬ 
ing optimal value of Problem (P2) by d^-^(A^ 2 )- Then, we 
find the optimal N 2 as = argmax G^min(^ 2 ) via 

one-dimensional search over N 2 . 

Eirst, we show one property of N 2 in the following propo¬ 
sition, which helps reduce its search space. 

Proposition 1: Eor Problem (P2), it holds that A"! < 
i.e., there should be no more points allocated to the inner ring 
than to the outer ring. 

Proof: Please refer to Appendix ■ 

With Proposition the set of potentially optimal values 
of N 2 is reduced to{l,2,...,y}, and the complexity of the 
one-dimensional search for N 2 is thus reduced by half. 

Then, let 002 {N 2 ) and P 2 (^ 2 ) denote the solution to Prob¬ 
lem (P2) with a given N 2 G {1,2,...,^}- Notice that 
with given N 2 , co’|(A' 2 ) can be first obtained by solving the 
following problem: 

(P2.1) min Ci^ 2 {N 2 ,uj 2 ). 

0<u;2<27r ’ 

Problem (P2.1) can be optimally solved by Algorithm 1 shown 
in Table |I] at the top of next page, for which the details are 
given in Appendix 


^Note that for the case of L = 1, the solution of (PI) can be easily shown 
to be the conventional A/'-ary PSK. For consistency, we also consider this case 
as a two-ring APSK constellation. 

^In the sequel, we follow the convention to assume that N is an integral 
power of two, thus is an even integer. 
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TABLE I 

Algorithm 1: Algorithm for solving Problem (P2.1) 

Input: N, N2 

Output: 002(^2), Cl^2(^2) 

Set Ni=N - N2. Initialize A" = |-|^|. 
for m = 0 to A^i — 1 do 
for n = 0 to A^2 — 1 do 

if -It < ^ - ^ < 0 then 

I = 

end 

end 

end 

Set K = \X\. Sort the elements in A' in a descending 
order: X(i) > X( 2 ) 

Set k* = argmin 

l<k<K-l 

^*(iV2) = , 

Cl2{N2) = COS 


> ... > X^K)- 
2 

*)-2^(fc*+l) . 


Let 2 (^ 2 ) denote the optimal value of Problem (P2.1). 


Note that in Table I since X, 


0 " 




■V 


(fe+i) 


( 0 , 


e 0 ,# 


for \/k e 


{1,2 ,K — 1}, we have CJ 2 (^ 2 ) ^ cos (^ 7 ^^ 5 Next, 
we obtain dj^in(^ 2 ) by finding p 2 (^ 2 ) through solving the 
following problem: 


(P2.2) 


max mm < 

^<P2<1 


P2 Vm, ^l + pl-2p2Cl^{N2)] . 


Problem (P2.2) is still non-convex with respect to p 2 . 
To find its optimal solution, w e observe that the inter-ring 
MED + P 2 ~ ‘^p 2 Cl 2 {^ 2 ) first decreases with p 2 when 
p 2 < C'l 2 (^ 2 ), and then increases with p 2 when 2 (^ 2 ) < 
P 2 < 1- Furthermore, the intra-ring MED of the inner ring, 
p 2 \/2B2, strictly increases with p 2 . Based on the above results, 
we discuss the solution of Problem (P2.2) in the following two 
cases. 

• Case 1: <^ 12 (^ 2 ) < In this case, we have p 2 > ^ > 
Cl 2 (^ 2 ), and the objective function of Problem (P2.2) is a 
non-decreasing function of p 2 . Therefore, we have 


pI{N 2) = 1 , 

Ci„(iV2) = min { v^, ^2 - 2Cl^{N2)] . 

( 22 ) 

This is consistent with our intuition that as ^ becomes large, 
the two rings converge and allocating all the signal points 
on the outer ring is optimal. 

• Case 2: Cl 2 (^ 2 ) > In this case, we further divide the 
feasible range of p 2 into [^5 2 (A^ 2 )] and (CJ 2 (^ 2 ), l], 

which are referred to as region I and region II, respectively. 
We then investigate the locally optimal p 2 within the two 
regions, which are denoted by p^ i(A^ 2 ) and p^ ii(A^ 2 ), re¬ 
spectively. Let d^in,i(^ 2 ) and d^in, 11 (^ 2 ) denote the local 
optimum of Problem (P2.2) for the two regions, respectively. 


The globally optimal solution to Problem (P2.2) is then 
given by 


pI{N2) 


PUN 2 ) 

P 2 11 (^ 2 ) of her wise. 


(23) 


Furthermore, for region I, since yl + ~ 2 p 2 C'i 2 (^ 2 ) 

and p 2 \/ 2 B 2 are monotonically decreasing and increasing 
with p 2 , respectively, the shape of the objective function of 
Problem (P2.2) as well as the corresponding ^ 21 (^ 2 ) are 
dependent on the values of Bi, B 2 , Cl^ 2 {^ 2 ) and We 
thus propose Algorithm 2 in Table |n to find ^ 2,1 (^ 2 ). for 
which the details are given in Appemix 


TABLE II 

Algorithm 2: Algorithm for finding p * i(A^2) 



On the other hand, for region II, similar to Case 1 above, 
the objective function of Problem (P2.2) is a non-decreasing 
function of p 2 since p 2 > 2 (^ 2 )- Therefore, we have 

P*2AN2) = 1, 

d*min,11(^2) = min { v^, V^2, ^2-2CI^{N2)] . 

(24) 

So far, we have found the optimal solution for ujI{N 2 ), 
pI{N 2 ) and the corresponding maximum MED d^in(^ 2 ) 
for any given N 2 . With Proposition the optimal number 
of points on the inner ring is then obtained as = 
argmax dj^in(^ 2 )- This completes our proposed algo- 
Ar2e{i,2,...,f } 

rithm for Problem (P2), which is summarized as Algorithm 
3 shown in Table |I^ at the top of next page. For any given N 
and ^, this algorithm can be shown to find the optimal solution 
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to Problem (P2) with worst-case complexity of 0{N^ \ogN). 


TABLE III 

Algorithm 3: Algorithm for solving Problem (P2) 


Input: N, ^ 

Output: No, cco, po 
for Ar 2 = 1 to f do 

Compute Bi and B 2 . 

Obtain cc 2 (^ 2 ) and 2 (^ 2 ) by Algorithm 1. 
if Cl* 2 (^ 2 ) < i then ’ 

I Obtain pl{N 2 ) and d^in(^ 2 ) by 

else 

Obtain ^ 2,1 (^ 2 ) and d^in,i(^ 2 ) by Algorithm 2. 
Obtain piji{N 2 ) and d^in, 11 (^ 2 ) by (p(. 

I Pim = Phm-, = Ci„,i(A'2). 

else 

P*2{N2) = pluiN2y, 

d*minm=d* 

end 

end 

end 

Set N 2 = arg max ^ 

Ar2e{i,2,...,f } 

^2*=^2*(A^2*)- 


^min,11(^2)- 


,(7V2);p*=p*(7V*; 


To illustrate our proposed design for A^-ary APSK constel¬ 
lation with 1/ = 2, we take the example of = 16 and 
show the optimal p 2 and A^| versus ^ in Fig. In addition, 
we provide in Fig. the constellation diagrams of optimal 
16-ary APSK with L = 2, for the cases of = 0.4, 0.5 and 
0.7, respectively. Furthermore, the optimal set of A^-ary APSK 
constellations with L = 2 as well as their maximum MEDs 
(denoted by for the case of A^ = 16 are shown in 

Table W at the top of next page, while those for the cases of 
A^ = 8, 32 and 64 are provided in Tables V-VII in Appendix 
D. It can be observed that for each N, the values of ^ are 
divided into several regions, each of which is associated with 
an optimal constellation design]^ We also show in Tables IV- 
VII the probability of each region for the case of i.i.d. Rayleigh 
fading channel with M = 2 or M = 4 transmit antennas, 
which is obtained through 10^ Monte Carlo trials. Moreover, 
for the cases of A^ = 2 and A^ = 4, it can be shown that for any 
^ G [0,1], the optimal solutions to Problem (P2) correspond to 
BPSK (binary phase shift keying) and QPSK (quaternary phase 
shift keying) constellations, respectively, with 2 — 

^min 4 = respectively. For other values of N, the optimal 
constellation designs can be obtained similarly. 


V. SuBOPTiMAL Solution to Problem (P2) 

Clearly, the amount of memory required to store the set 
of optimal A^-ary two-ring APSK constellations at the trans¬ 
mitter/receiver increases with the number of regions resulting 

^Particularly, the optimal constellation parameters as well as le 
constant over Regions 1, 3, 5 and 7. This is because for each of these regions, 
the optimal constellation corresponding to the left boundary value of ^ is 
also feasible for the other values of ^ in the same region. 



R 


Fig. 3. P 2 and versus ^ when = 16. 



Fig. 4. Constellation diagrams of optimal 16-ary APSK with L = 2 for 
different 


from the optimization, each corresponding to a different con¬ 
stellation. Moreover, the probability of each region used in 
fading channel may also vary significantly from one another, 
e.g., with A^ = 16 and i.i.d. Rayleigh fading channel, it can 
be observed from Table |IV] that the probability of Region 
1 is much higher than that of Region 2 for the cases with 
M = 2 and M = 4 transmit antennas, respectively. This 
thus motivates us to design a suboptimal constellation set for 
N-Sivy two-ring APSK with a smaller number of regions or 
constellations. From Tables IV-VII, we have the following 
observations: 

• First, for each N, the constellation corresponding to Region 
1 has the maximum MED among all the constellations; 
moreover. Region 1 generally has the largest probability, 
and the probability increases as N increases; 

• Second, for each N, the constellation corresponding to the 
last region is always feasible for the values of ^ in all the 
other regions. 

Therefore, we propose a suboptimal constellation set design 
with no more than two regions. Specifically, in the suboptimal 
design, the range of ^ and the constellation corresponding to 
the first region are identical to those in the optimal set; while 
the second region (if any) covers the remaining range of ^ 
and uses the constellation for the last region in the optimal 
set|^ As an example, in the proposed suboptimal design for 
N = 16, the constellation for ^ G [0, 0.4603) is the same as 
that for Region 1 in Table while the remaining range of 


^For cases such as N = 2 and N = A where the optimal set has a single 
region, the suboptimal set is identical to the optimal set. 
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TABLE IV 

Optimal constellation set eor 16- ary APSK with L = 2 


Region 

r 

R 

P2 

iV| 

^2 

^min,16 

Prob^=^ 

Prob^'^=^ 

1 

[0,0.4603) 

0.4603 

5 

0.01827r 

0.5411 

0.7596 

0.9994 

2 

[0.4603,0.4839) 

r 

R 

5 

0.01827r 

iaf -2{i)CU(Ni) + l 

0.0245 

0.0002 

3 

[0.4839,0.5176) 

0.5176 

4 

0.08337r 

0.5176 

0.0323 

0.0002 

4 

[0.5176,0.5588) 

r 

R 

4 

0.08337r 

\/(K)'-2(ac;,,(wj) + i 

0.0352 

0.0001 

5 

[0.5588,0.6302) 

0.6302 

8 

0.12507r 

0.4824 

0.0505 

0.0001 

6 

[0.6302,0.8477) 

r 

R 

8 

0.12507r 


0.0844 

0.0000 

7 

[0.8477,1] 

1 

8 

0.12507r 

0.3902 

0.0135 

0.0000 


^ G [0.4603,1] uses the constellation for Region 7 in Table 

lYl 

VI. Numerical Results 

In this section, we provide numerical results. Consider a 
single-user MISO channel modeled by , where 

represents Rayleigh fading channel coefficients with i.i.d. 
elements hi ^ CN'{0, 1), and p = —90dB denotes the channel 
power attenuation due to path lossThe average receiver noise 
power is set to be = — 94dBm[^ The average signal-to- 
noise ratio (SNR) is thus defined as SNR = ^. All the results 
below are averaged over 10^ independent channel realizations. 

A. Fixed-Rate Adaptive Receiver Constellation 

In this subsection, we evaluate the performance of CE 
precoding with fixed-rate adaptive receiver constellation, based 
on the optimal set of two-ring APSK constellations or the 
suboptimal alternative. Moreover, we examine the effect of 
imperfect CSIT on the performance of the proposed schemes. 

1) Performance of Optimal Constellation Set for IQ-ary 
Two-Ring APSK: Consider the case of = 16, and we adopt 
the optimal constellation set for 16-ary APSK with L = 2 (as 
shown in Table |IV]). In addition, we consider two benchmark 
schemes for performance comparison: 

• Benchmark Scheme 1: CE precoding with fixed rectangular 
16-QAM constellation; 

• Benchmark Scheme 2: CE precoding with adaptive receiver 
constellation, by using rectangular 16-QAM constellation 
when it is feasible, i.e., ^ and 16-PSK constellation 
otherwise. 

In Eig. we show the average SER comparison of our 
proposed scheme versus the two benchmark schemes for the 
case of M = 2 or M = 4 transmit antennas, respectively. 

Eor both setups, it is observed from Eig. that Benchmark 
Scheme 1 results in severe error floor due to its infeasibility 
when ^ > 1 It is also observed that our proposed scheme 

^We consider a path loss model given by /3 = K |l^ , where we 

assume K = 10“^, a = 3, do = Im and d = 100m. 

^^We assume the bandwidth of the transmitted signal is lOMHz, the power 
spectral density of the receiver noise is — 174dBm/Hz, and the noise figure 
due to receiver processing is lOdB. 



(a) M = 2 



Fig. 5. Average SER comparison of fixed-rate CE transmission schemes. 

outperforms Benchmark Scheme 2 by 1.42dB and 1.52dB in 
SNR at the given SER of 10“^ for the cases of M = 2 
and M = 4, respectively]^ This performance gain is due to 
our optimization of constellation tailored for CE transmission: 
when ^ <1 , the MED of our proposed optimal constellation 
is 0.5411, while that of rectangular 16-QAM is 0.4714; when 

Since rectangular 16-QAM has larger MED than 16-PSK (0.4714 versus 
0.3902), it can be further shown that our proposed scheme outperforms CE 
precoding with fixed 16-PSK constellation at the receiver, which is always 
feasible regardless of the channel realization. 


































^ the MED of 16-PSK is also no larger than that of our 
proposed optimal constellations. 

2) Performance of Suboptimal Constellation Design: In 
Fig. we consider M = 2 and show the average SER com¬ 
parison of CE precoding with adaptive A^-ary two-ring APSK 
receiver constellation based on the optimal and suboptimal 
constellation designs, for the case of = 16 or A^ = 64. 
It is observed that for both constellation sizes, the SNR gap 
between the optimal and suboptimal designs at any SER is 
very small and almost negligible. Furthermore, this SNR gap 
is observed to decrease as N increases, since the probability 
of Region 1 increases with increased N. The above results 
indicate that the proposed suboptimal constellation design 
of low complexity is practically efficient for fixed-rate CE 
transmission. 



Optimal design; M=2 
Suboptimal design; M=2 

- EOT w/rectangular 16-QAM; M=2 
-Optimal design; M=4 

Suboptimal design; M=4 

- EOT w/ rectangular 16-QAM; M=4 
' Perfect CSIT; Optimal design; M=2 

Perfect CSIT; Suboptimal design; M=2 
Perfect CSIT; EOT w/ rectangular 16-QAM; M=2 

- Perfect CSIT; Optimal design; M=4 
Perfect CSIT; Suboptimal design; M=4 

- Perfect CSIT; EOT w/ rectangular 16-QAM; M=4 
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8 10 12 
SNR*' (dB) 


Fig. 7. Average SER versus average channel training SNR. 



Fig. 6. Average SER comparison of the optimal and suboptimal constellation 
designs for A^-ary APSK with L = 2. 

3) Effect of Imperfect CSIT: To realize CE precoding with 
the proposed fixed-rate adaptive receiver constellation, the 
instantaneous values of r and R are required to be known 
at the transmitter to select the constellation, which is possible 
since the exact channel coefficients are assumed to be 
known at the transmitter to compute the transmitted signal 
phases. In this part, a time-division duplex (TDD) system is 
considered, where the transmitter acquires CSI via the reverse 
link channel training by the receiver under the assumption of 
channel reciprocity!^ In practice, perfect channel estimation 
is generally unavailable at the transmitter, which motivates us 
to examine the effect of imperfect CSIT on the performance 
of our proposed schemes. 

Let h denote the minimum mean-square error (MMSE) 
estimate of the channel vector h at the transmitter, h 
is modeled by h = h — Ah, where Ah denotes the estimation 
error. Let SNR^^ denote the average SNR of the training signal 
received at the transmitter in the presence of i.i.d. CSCG noise. 
The distribution of each element in Ah can be shown to be 
given by Aft,j ~ CM{Q, i_|_gNRtr )i * = i.e., 

larger SNR^^ leads to more accurate CSI estimation. In Fig. 

the other hand, applying frequency-division duplex (EDD) to CE 
precoding schemes requires the forward link training signal at each transmit 
antenna to be designed under the CE constraint. This is a challenging problem 
that has not been addressed yet in the literature to the best of our knowledge, 
thus is left as a possible direction for future work. 


1^ we show the average SER versus SNR^^ for our proposed 
CE precoding with adaptive receiver constellation based on the 
optimal and suboptimal designs of 16-ary two-ring APSK for 
the case of M = 2 or M = 4 transmit antennas, respectively, 
where we assume N = 16 and SNR = 20dB for the purpose 
of exposition. For each scheme, we also show the average SER 
with perfect CSIT as benchmarks. Under both setups, it is ob¬ 
served that the SER can be effectively improved by increasing 
SNR^^, which shows that the accuracy of CSIT is crucial to 
the performance of our proposed schemes. This is because 
the nonlinear CE precoding and adaptive receiver constellation 
selection in our proposed scheme are both based on CSIT. In 
addition, we show in Fig. the average SER versus SNR^^ 
and that with perfect CSIT for equal gain transmission (EGT) 
with fixed rectangular 16-QAM constellation, which is a linear 
non-CE precoding scheme. It is observed that the performance 
of our proposed schemes is almost equally sensitive to CSIT 
inaccuracy as that of EGT with fixed rectangular 16-QAM. 

B. Variable-Rate Adaptive Receiver Constellation 

In this subsection, we consider the case of variable-rate 
transmission with CE precoding. In general, the aim of 
variable-rate adaptive receiver constellation design is to max¬ 
imize the spectral efficiency log 2 N (in bps/Hz) by jointly 
adapting N and S to the instantaneous CSI, such that the N- 
ary constellation S is feasible at the receiver, and the receiver 
SER is lower than a target value Pq. Similar to the fixed-rate 
CE transmission case, for simplicity, we replace the constraint 
on the exact receiver SER with that on the union bound of 
SER. Furthermore, for any given ^, we assume S to be within 
the family of optimal A-ary two-ring APSK constellations 
with N e ff, where ff denotes the set of possible values of N. 
Therefore, we formulate the following optimization problem 
with given r, R, and M for the case of variable-rate CE 
transmission for average spectral efficiency maximization: 

(P3) max logo N 

s.t. {N - \)Q < Pe, (25) 
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where ^ can be readily obtained by solving Problem (P2) 
with Algorithm 3 for the given The optimal solution to 
Problem (P3) can be obtained via one-dimensional search over 

in a decreasing order, namely, from N = maxA^ to = 

_ ' NeAT 

min N, until (25) is satisfied. 

NeAf 

We consider me case of A/" = {2,4,8,16,32,64} and Pe = 
10“^ to evaluate the performance of our proposed variable- 
rate CE transmission scheme. For performance comparison, 
we also consider the following benchmark scheme: 

• Variable-rate transmission with CE precoding based on the 
family of rectangular QAM constellations: the constellation 
size for given r and R is chosen as the maximum N G 
JV such that the A^-ary rectangular QAM is feasible at the 
receiver and yields an SER union bound below Pep^ 



(a) Proposed scheme with two-ring APSK constellation 



(b) Benchmark scheme with rectangular QAM constellation 
Fig. 8. Optimal constellation size N* versus ^ and R. 

For illustration, the optimal constellation size A/"* versus 
^ and R for both the proposed and benchmark schemes are 
provided in Fig.[^ Note that we let N* = 1 (i.e., log 2 A/"* = 0 
in the objective function of (P3)) denote the case of no 
transmission, i.e., the SER constraint in ( [25] ) cannot be met 
for any N e JV. In Fig. we compare the average spectral 

^^There are various constellation designs for rectangular 8-QAM with 
different MEDs and thus different expressions for the SER union bound. To 
avoid ambiguity, we use the one shown in Figure 4.3-7 (c) in 0. 


efficiency of our proposed scheme versus the benchmark 
scheme, for the case of M = 2 or M = 4, respectively. 



(a) M = 2 



Fig. 9. Average spectral efficiency comparison of variable-rate CE transmis¬ 
sion schemes. 

For the case of M = 2, it is observed from Fig. (a) that 
our proposed scheme outperforms the benchmark scheme by 
1.56dB in SNR at the average spectral efficiency of 3bps/Hz. 
It is also observed that the SNR gap between the two schemes 
increases as the desired average spectral efficiency grows. By 
examining the feasibilities and MEDs of rectangular QAM 
constellations, this performance gain can be explained as 
follows. First, for N G {2,4,8,16}, the MED of the A^-ary 
rectangular QAM constellation can be shown to be no larger 
than that of the optimal A^-ary two-ring APSK constellation. 
Furthermore, for N G {8,16,32,64}, the probability of the 
A^-ary rectangular QAM constellation being infeasible can be 
shown to be nonzero and increases as N grows. 

For the case of M = 4, it is observed from Fig. (b) that 
the SNR gain of our proposed scheme over the benchmark 
scheme is reduced to 1.39dB at the average spectral efficiency 
of 3bps/Hz. Furthermore, the benchmark scheme even outper¬ 
forms our proposed scheme by 0.72dB of SNR at the average 
spectral efficiency of 5bps/Hz. This is due to the following 
reasons. First, for N G {8,16,32,64}, the probability of 
the N-nry rectangular QAM being infeasible is much lower 
for the case of M = 4 than that for the case of M = 2, 
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e.g., rectangular 16-QAM is infeasible for 0.3% of channel 
realizations when M = 4 instead of 40% when M = 2 
under the assumed i.i.d. Rayleigh fading channels. Second, 
when rectangular 32-QAM and 64-QAM constellations are 
feasible, their MEDs (0.3430 and 0.2020, respectively) are 
larger than those of the optimal 32-ary and 64-ary two-ring 
APSK constellations (0.3238 and 0.1778, respectively). Note 
that the smaller MEDs of our proposed APSK constellations 
are due to the limited number of concentric rings (i.e., L = 2), 
which can be improved by finding the solution to Problem 
(PI) with L > 2, especially for large values of N. This is an 
interesting problem to be investigated in our future workp] 


Assume S and S' have the same inner ring radius P 2 (^ 2 )- 
Eor I G {1,2}, let d/ and dj denote the intra-ring MEDs of 
the Ith ring for S and S', respectively. Let di ^2 and d} 2 denote 
the inter-ring MEDs for S and S', respectively. The MEDs of 
S and S' are thus given by dmin = min{(ii,^ 2 ,^ 1 ^ 2 } and 
^min = 2 }^ respectively. Erom (13), we find 

that d 2 < p*^n 2 ) ~ < p 2 {^ 2 )di = d' 2 . Moreover, it 

can be observed from (14) that di^2 = 2 - Therefore, we have 

dmin < ^2 < mm{d]Jd' 2 } and dmin < g^i ,2 = <^ 1 , 2 ’ which 
yield d^iin ^ ^min’ that d^iin — ^min('^ 2 ) and ^ 
^min(^-^ 2 ), we thus have d^-^{N 2 ) < d^-^{N-N 2 ). This 
completes our proof of Proposition 


VII. Conclusions 

This paper studied the design of adaptive receiver constel¬ 
lation for CE precoding in a single-user MISO fiat-fading 
channel with arbitrary number of transmit antennas. Eor the 
case of fixed-rate transmission, an efficient algorithm was 
proposed to obtain the optimal APSK constellation with two 
rings that is both feasible and of the maximum MED for any 
given constellation size and instantaneous CSI. It was shown 
that our proposed scheme significantly improves the SER 
performance of CE precoding with fixed receiver constellation 
in fading channel. In addition, a suboptimal design for fixed- 
rate adaptive two-ring APSK constellations was proposed, 
which is of low complexity and yet achieves near-optimal 
SER performance. Eurthermore, we proposed a variable-rate 
CE transmission scheme based on the family of optimal fixed- 
rate adaptive two-ring APSK constellation sets and examined 
its average spectral efficiency as compared to that based on 
the rectangular QAM constellations. 

An appealing direction of future work is to study the optimal 
adaptive APSK constellation design with an arbitrary number 
of concentric rings, which is expected to further enhance the 
performance of the proposed design with two rings. It is also 
interesting to extend the results in this paper to the more 
general setups with multiple users/receive antennas. 

Appendix A 

Proof of Proposition[T] 

We prove Proposition Qby showing that for any A ^2 > y, 
it is true that d^in(^ 2 ) < ^^min(^ - ^ 2 ). Let cc|(A^ 2 ) and 
p 2 {N 2 ) denote the solution to Problem (P2) for any given 
N 2 G {1,...,A^ — 1}. Consider two sets of A^-ary two- 
ring APSK constellations denoted by S and S', respectively. 
Assume the numbers of points on the inner ring of S and 
S' are given by A ^2 ^ {t + ..., — l} and N — N 2 , 

respectively. Assume the phase offsets of the inner ring of 
S and S' are given by cc5(A^2) and 27r — (A^ 2 ), respectively. 

^"^It is worth noting that there are two new challenges for solving Problem 
(PI) with L > 2 compared to the case with L = 2. First, note that the terms 
with \l — h\ / 1 denote the inter-ring MEDs between any two signal 
points that are located on two non-adjacent rings, which will be present in 
the case of L > 2. Additional effort is thus needed to deal with the constraint 
in (7) for the case of L > 2. Second, it can be shown that the feasible set of 
{^l}l ^=2 of L > 2 has cardinality of Therefore, finding 

the optimal with affordable complexity is also a difficult problem 

for large values of N and L. 


Appendix B 

Algorithm for solving Problem (P2.1) 

We first observe that by symmetry, the constraint of Problem 
(P2.1) can be modified as 0 < u ;2 < which yields -b 

W 2 - ^ e j^-2Tr,2TT-j^ + j^ C {-2tt, 2-k]. Note 
that for (p G (—27r,27r], coscj) monotonically increases with 
1101 — 7r|. Therefore, Problem (P2.1) is equivalent to 


min max f{m^n^uj 2 ), 

0<cu2<;^ rnGXi,nGX2 


(26) 


where Ii = {0,1,..., Vi — 1}, I 2 = {0, l!”-i-V 2 — 1} and 


f{m,n,uj2) = 


0 J 2 


( 27rn 27rm A 

N 2 Ni ) 


Eor the case of A^ > 2, we have the following results for 
f{m,n,uj 2 ): 

f{m,n,uj2) = 

27rn _ 277771 ^ r —277 qI 

N 2 Vi ^ L Ni 

< max{/(l, 0, 0 J 2 ), /(0,0, W 2 )}, 


TT + W2 + ^ - 


LO 2 ■ 


27777 277771 

N 2 Ni 

( 27777 277771 A 

N 2 Ni ) 


27777 _ 277771 

Ni 


, , I Z77n 

‘^2 + tv2-- 


— TT 


27777 

N 2 


277771 PI 
Ni ^ 

(27) 

as 


f 27rn 

27rm^ 


Vi ) 


Based on (27), we further simplify Problem 


0<^2<^ 771GXi,77GX2 

^ 27^ _ 27rm ^ f —27r <-.1 

N2 ^^1 ^ [ Ni 

(28) 

To solve Problem ([28]), we first select the set A’ = 

and sort the el- 


ements therein as > X( 2 ) > ••• 


> 


(A), 


where 


K = |A'|. Next, we find the locally optimal solution when 
CC 2 G [—for each/c G {1, 2,..., iZ —1}, which 
is given by 002 (^ 2 , k) = Einally, the globally 

optimal solution to Problem (28) is given by (^ 2 (^ 2 ) = 




Ufc*) where /c* = argmin 

l<k<K-l 


-X, 


(fc) 


. The cor¬ 


responding optimal value to Problem (P2.1) is thus given by 

. Eor the case of A^ = 2, i.e.. 




(fc*+i) 


Cl* 2 (^ 2 ) = cos ( 

A^i = A ^2 = it can be easily shown that uj 2 {N 2 ) = and 
C{ 2 (^ 2 ) = cos consistency, we include this case 

in Algorithm 1 by initializing A’ = An illustration 
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of Algorithm 1 is shown in Fig. 10 for the case of A" = 16 
and N 2 = 5. 


Ad ^(2) "^(3) ^(4) ^(5) ^(6) 



Fig. 10. Illustration of Algorithm 1 with N = 16, N 2 = 5. 


Appendix C 

Algorithm for finding 1 (^ 2 ) 

We find ^ 21 (^ 2 ) by examining the objective function 
of Problem (P2.2), which is characterized by the inter- 
se ction points of the lines given by \/2Bi, p 2 \/ 2 B 2 , 
^1 + ^ 2 “ ‘2.Cl 2 {^ 2 )p 2 as well as the interval boundaries 
given by p 2 = ;^ and p 2 = 2 (^ 2 ). First, note that if 

Ni < 2 , we have ^^ 2 (^ 2 ) < 0 < which contradicts 
the assumption of Case 2 in Section IV. Thus, we only need 
to consider Ai > 3 in the sequel, which yields 2 (^ 2 ) > 

cos (^ 7 ^^ ^ Based on this result, it can be proved that 
when P 2 = Q, 2 (-^ 2 ), \J^+p2- 2 Cr. 2 (-^ 2)/32 < < 

1+pI-2CI^{N2)p2 < < \/2Br. 


p 2 y/ 2 B 2 and 

T hen, we denote the (lef t) intersection point of p2\/2B2 and 
+ pI- 2 Ci* 2 (^ 2 )P 2 by (P 2 , P 2 V 2 B^), where p 2 is given 

by 


P2 — 


C'r,2(^2)-VC'r,2(^2)+2H2-l 
1 - 2 S 2 

1 


2C* 


^1,2(^2) 

Moreover, we find 
between \/ 2 Bi and 


B2 = \. 


(29) 


the 


_ (left) intersect ion point 

^1 + pI- 2C^ 2 i^ 2 )p 2 given by 

(Ci*2 {N2)- yJCt%{N2)- 1 + 2 Bi, V2S7), and the 
intersection point between \/2Bi and P 2 V 2 B 2 given by 
(\/^’ Finally, based on the above results, we 

discuss the solution of ^ 21 (^ 2 ) in fonf cases, which are 
summarized in Algorithm 2 shown in Table Note that if 
there are multiple optimal solutions, we choose the one to 
maximize the SNR. For the purpose of illustration, examples 


of the four cases are shown in Fig. 11 for the case of A = 16. 



(a) Case i, with N 2 = Q and ^ =0.25 



(b) Case ii, with N 2 = 4: and ^ =0.25 



- 

- P2V2B2 


- - ^l+pl-2p2Cl2{N2) 

' - ^ {p2,P2\/2.B'^ 

y 

0.DD_^^ ^ 

D--.... 

:(p;.l(A^2),<in.l(iV2)l 

C{p{N2) 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

P2 


(c) Case iii, with N 2 = 4 and -£ = 0.55 



(d) Case iv, with N 2 = 4: and ^ = 0.45 


Fig. 11. Illustration of Algorithm 2 with N = 16. 


Appendix D 

Optimal Constellation Sets for 8,32 and 64-ary 
APSK WITH L = 2 
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TABLE V 

Optimal constellation set eor 8-ary APSK with L = 2 


Region 

r 

R 

P2 



^min,8 

Prob^=‘^ 

Prob^=^ 

1 

[0,0.1495) 

0.1495 

1 

0.14297r 

0.8678 

0.2925 

0.9766 

2 

[0.1495,0.2705) 

r 

R 

1 

0.14297r 


0.2117 

0.0170 

3 

[0.2705,1] 

1 

4 

0.25007r 

0.7654 

0.4958 

0.0064 


TABLE VI 

Optimal constellation set eor 32- ary APSK with L = 2 


Region 

r 

R 

P2 


^2 

^min,32 

Prob^=^ 

Prob^=‘‘ 

1 

[0,0.6764) 

0.6764 

13 

0.093l7r 

0.3238 

0.9282 

9.9999 X 10-1 

2 

[0.6764,0.6873) 

r 

R 

13 

0.093l7r 

\/(i)'-2G*2(N|)i + l 

0.0054 

2.5000 X 10-® 

3 

[0.6873,0.6902) 

0.6902 

12 

0.01677r 

0.3129 

0.0014 

0.7000 X 10-^ 

4 

[0.6902,0.7074) 

r 

R 

12 

0.01677r 

\/(K)“-2C,-,2(iV2')s + l 

0.0079 

3.0500 X 10-® 

5 

[0.7074,0.7583) 

0.7583 

16 

0.06257r 

0.2959 

0.0200 

4.2600 X 10“*^ 

6 

[0.7583,0.9615) 

r 

R 

16 

0.06257r 

\/(i)"-2CT.2(«{)i + i 

0.0363 

1.8700 X 10“® 

7 

[0.9615,1] 

1 

16 

0.06257r 

0.1960 

0.0008 

0.0000 


TABLE VII 

Optimal constellation set eor 64- ary APSK with L = 2 


Region 

r 

R 

P2 


^2 

^min,64 


ProL'i=^ 

1 

[0,0.8222) 

0.8222 

29 

0.00307r 

0.1778 

0.9811 

999.9998 x 10“® 

2 

[0.8222,0.8257) 

r 

R 

29 

0.00307r 

\/(s)‘‘-2Cr,2(W2-)B + l 

0.0008 

0.0000 

3 

[0.8257,0.8261) 

0.8261 

28 

o.ouqtt 

0.1743 

0.0001 

0.0000 

4 

[0.8261,0.8321) 

r 

R 

28 

o.ouqtt 


0.0013 

0.0300 X 10-® 

5 

[0.8321,0.8584) 

0.8584 

32 

0.03127r 

0.1683 

0.0051 

0.1000 X 10“® 

6 

[0.8584,0.9903) 

r 

R 

32 

0.03127r 

V(i)'-2G*2(N|)i + l 

0.0115 

0.0300 X 10-® 

7 

[0.9903,1] 

1 

32 

0.03127r 

0.0981 

0.0001 

0.0000 
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